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Strong Coupling Model for String Breaking on the Lattice. 
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A model for SU(n) string breaking on the lattice is formulated using strong coupling ideas. It gives an explicit 
picture of string breaking, in the presence of dynamical quarks, as a mixing process between a string state and 
a two-meson state. An analysis of the Wilson loop shows that the evolution of the mixing angle as a function 
of separation may obscure the expected crossover effect. If a sufficiently extensive mixing region exists then an 
appropriate combination of transition amplitudes may help to reveal the effect. The sensitivity of the mixing 
region to the values of the meson energy and the dynamical quark mass is explored. 
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1. Introduction 

Much effort is being devoted to dynamical 
quarks [[Ir[|. For a review see S Giisken §. For 
currently feasible values of the quark mass it is 
not expected that there will be dramatic effects 
on the computed hadron spectrum. String break- 
ing is an effect that is only possible in the pres- 
ence of dynamical quarks. Here we show explic- 
itly that string breaking does occur on the lattice 
and shows up as a mixing phenomenon. 

The physical situation envisaged in lattice cal- 
culations comprises a static quark and anti-quark 
separated by a spatial distance of R lattice units. 
The two static particles may either support a 
gauge string that runs between them or sepa- 
rately bind a quark and an anti-quark to create 
a two-meson state. When the string is stretched 
sufficiently, its energy coincides with the energy 
of the two meson state. In the neighbourhood of 
this critical separation the two different physical 
states mix permitting the occurence of transitions 
between them. 

2. Strong Coupling Model 

We recall the rules for evaluating simple graphs 
in the strong coupling limit of SU (n) gauge the- 
oryil- 

1. A factor of (/3/2n 2 ) for each plaquette. 

2. A factor of 2n (1 + j.e) /2 for each Wilson 
quark line in the direction of the unit vector 
e . Here k is the standard quark hopping 



parameter. 

3. Factors of (—1) and 1 jn and a trace over the 
spin matrix factors for each internal quark 
loop. 

2.1. Simple String Model 

In the absence of dynamical quarks the corre- 
lation function of a string of length R over an 
imaginary time interval T is the standard R x T 
Wilson loop. In leading strong coupling approxi- 
mation the above rules give for this string-string 
propagator 



Gss(T) 



-am 



(1) 



where the dimensionless string tension a = 
— log(/3/2n 2 ) . The energy of the string state 
is V{R) = aR . 

2.2. Model for Mesons 

Our model for mesons is one in which a light 
quark propagates along a static (anti-)quark line 
with a hopping parameter k! that encodes the en- 
ergy of the static bound state. Using this rule, we 
find that the meson propagator has the structure 
g(T) = (1 + 7o)/2 (2n') T . In other parts of the 
diagrams of the model the hopping parameter re- 
tains the value k related to the light quark mass. 

The propagator for two mesons moving inde- 
pendently, each bound to its static quark, is 



gW(T)®g^(T) 



'7o 



(i) 
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In fact only a particular combination of quark- 
anti-quark spin wave functions is relevant to the 
mixing phenomenon. It is obtained by completing 
the light quark loop with the matrices (l±7i)/2 , 
and including a factor of (—1) . The two-meson 
propagator becomes 



'mm (T) = \ ((2k') 2 ) T 



(3) 



We have then (2k') 2 = e Em , where Em is the 
combined energy of the two bound mesons. 

3. String Breaking 
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Figure 1 . Wilson loop (heavy line) containing in- 
ternal quark loops (light lines). 



Our model for string breaking involves sum- 
ming over planar graphs that incorporate tran- 
sitions between the string and meson states. A 
typical graph is shown in Fig 1. The correspond- 
ing contribution to the Green's function is 



[e-^] Tl (2k') 



-am T 3 



t 
n 

where 
t = Tr 



(2«) 



7o 



W] 4 



l + 7i 



(2k') 2 



(2k) 



7o 



^-aRl Tl 
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The factors involving k (as opposed to k') are as- 
sociated with the vertical sides of the quark loops 
and the factor V\V, where W is an energy, mea- 
sures the rate of transition from a meson state to 
qQ pair. (This factor was omitted in the original 
formulation |To[| .) Note that the trace of the in- 
ternal quark loop, which has the value t = — h, 
is equivalent to the quark-anti-quark spin projec- 
tion mentioned previously. 

The structure of our diagrams is as follows: 



1. A factor 



that propagates the 



string by one time step. 

2. A factor b = (2k') 2 that propagates the two- 
meson state by one time step. 

3. A factor c = (W/\f2n) (2k) R associated 
with the transition from string to two- 
meson state and vice- versa. 

To describe the transition from an initial time 
to time T we need a 2 x 2 matrix of transition 
amplitudes 

r(T\ - ( G ss{T) G S m(T) \ , s 

1 j ~ v. G MS (T) G M m(T) J ■ W 

then the above stepping procedure can be repre- 
sented by G(T + 1) = AG(T) , where the matrix 
A is given by 



A = 



be b 



(5) 



If for defmiteness we set G(0) = 1, then G{T) 
(A) T . We can set A = DOKO^D- 1 where 

a \ „ / cos 9 — sin 



D 



A : 



Vb 

X+ 
A_ 



,o 



suit 



cos( 



,(6) 
(7) 



The entries in A are the eigenvalues of A 

A± = i {(a + 6) ± y/(a - b) 2 + 4a6c 2 } . (8) 

The corresponding eigenvectors are the columns 
of DO. The mixing angle 9 is 



taS= -M+ J pw ^ (g) 

2v ab c 
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4. Transition Amplitudes 

We assume that each transition amplitude be- 
gins and ends with an appropriate propagator. 
We hnd for the Wilson loop 



Gss(T) = (1,0)A 



T-l 



(10) 



with the result 
Gss(T) = a (cos 2 9 AT" 1 



sin 2 6 A T_1 ) 



■ (11) 

Qmm{T) — b (sin 2 9 A^ _1 + cos 2 9 A 7 ^ 1 ) , (12) 



Similar reasoning yields 



Gms{T) = V^sinflcos^A^T 1 - A T ~ [ 1 
All amplitudes see both exponentials, e 



• (13) 

-E±T _ 

A^ . The the lower energy exponential should 
dominate the asymptotic behaviour. However 
what is observed will be influenced by the be- 
haviour of the mixing angle. When aR « Em 
then a » b and 9 ~ ,E + ~ aR ,£L ~ E M . 
Under these circumstances the coupling to the 
state with energy E_ will vanish and only the 
exponential associated with E + will be observed 
as expected. When aR >> Em then b >> a 
and 9 ~ § ,£L ~ B M , E_ ~ eri? . Now the 
coupling to the state with energy E + will van- 
ish and only the exponential associated with E_ 
will be observed. That is, both above and be- 
low the crossover point, only the original string 
behaviour e~ v ^ T , will be observed. The move- 
ment of the mixing angle therefore obscures the 
crossover phenomenon in which the string energy 
is expected to be bounded by the energy of the 
two-meson state. Complementary results hold for 
Gmm(T) . It is clear that Gsm(T) is suppressed 
outside the mixing region, where sin 9 cos 9 ~ . 

5. Mixing Region 

To see mixing effects directly in the Wilson loop 
the mixing region in R, for which both sin (9 and 
cos# are reasonably large, must be of sufficient 
size to be resolved on the lattice. The range of 
the mixing region in R can be estimated as 

/ it dR\ [2W m w ,„ / \ 

Ai ?= hr-rH =7r\ -—e m " EM/a .(14 
\2 d9 V 3 a y ' 



If we estimate the parameters in the model as 
a' 1 = 1.5 GeV, a 1/2 = 0.427 GeV, W ~ E M = 
1.0 GeV then the critical distance is R c = 8.22 
lattice spacings. For m q —.l GeV we then find 
Ai2=18.3 lattice spacings, while for m q =.5 GeV 
we obtain Ai?=2.0 lattice spacings. We see then 
that if the dynamical quark is light enough and 
the meson energy not too high there is a rea- 
sonable chance of observing the string breaking 
within the mixing region. The exponential fac- 
tor in eq(|l4|) guarantees that if the meson energy 
is larger than 1 GeV the mixing region will be 
rapidly restricted as the value of m q is raised. 
An appropriate combination of transition ampli- 
tudes may also help in identifying the mixing pro- 
cess. For example in the model the combination 
bGss{T)+aGMM(T) removes the influence of the 
mixing angle variation. 
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